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Abstract. We address the nonGaussianity (nG) of states obtained by weakly 
perturbing a Gaussian state and investigate the relationships with quantum 
estimation. For classical perturbations, i.e. perturbations to eigenvalues, we found 
that nG of the perturbed state may be written as the quantum Fisher information 
(QFI) distance minus a term depending on the infinitesimal energy change, i.e. 
it provides a lower bound to statistical distinguishability. Upon moving on 
isoenergetic surfaces in a neighbourhood of a Gaussian state, nG thus coincides 
with a proper distance in the Hilbert space and exactly quantifies the statistical 
distinguishability of the perturbations. On the other hand, for perturbations 
leaving the covariance matrix unperturbed we show that nG provides an upper 
bound to the QFI. Our results show that the geometry of nonGaussian states 
in the neighbourhood of a Gaussian state is definitely not trivial and cannot be 
subsumed by a differential structure. Nevertheless, the analysis of perturbations 
to a Gaussian state reveals that nG may be a resource for quantum estimation. 
The nG of specific families of perturbed Gaussian states is analyzed in some details 
with the aim of finding the maximally non Gaussian state obtainable from a given 
Gaussian one. 



1. Introduction 

NonGaussianity (nG) is a resource for the implementation of continuous variable 
quantum information in bosonic systems [T] . Several schemes to generate nonGaussian 
states from Gaussian ones have been proposed, either based on nonlinear interactions 

or on conditional measurements [3 II [SI II [3 [i [HI [ini [IS [II [II [II [111 113 III III 
[2n [TBI [22l [23l [24l [25] . In many cases the effective nonlinearity is small, and so it is the 
resulting nG. It is thus of interest to investigate the nG of states in the neighbourhood 
of a Gaussian state, i.e. the nG of slightly perturbed Gaussian states. Besides the 
fundamental interest [26] this also provides a way to assess different deGaussification 
mechanisms, as well as nG itself as a resource for quantum estimation. Indeed, in an 
estimation problem where the variation of a parameter affects the Gaussian character 
of the involved states one may expect the amount of nG to play a role in determining 
the estimation precision. 
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Quantum estimation deals with situations where one tries to infer the value 
of a parameter A by measuring a different quantity X, which is somehow related 
to A. This often happens in quantum mechanics and quantum information where 
many quantities of interest, e.g. entanglement |27[ 128] . do not correspond to a 
proper observable and should be estimated from the measurement of one or more 
observable quantities |29| . Given a set {q\\ of quantum states parametrized by 
the value of the quantity of interest, an estimator A for A is a real function of the 
outcomes of the measurements performed on g\. The quantum Cramer-Rao theorem 
[30l [3T| [32I [33] establishes a lower bound for the variance Var(A) of any unbiased 
estimator, i.e. for the estimation precision, Var(A) > {MH{X))^^ in terms of the 
number of measurements M and the so-called quantum Fisher information (QFI), 
which captures the statistical distinguishability of the states within the set. Indeed, the 
QFI distance itself is proportional to the Bures distance ^^(ei, Q2) = 2[1 — F{gi, ^2)], 
-^"(^17^2) = Tr [ y/ ^/q\Q2 ] being the fidelity, between states corresponding to 
infinitesimally close values of the parameter, i.e., in terms of metrics, i?(A) — ^g\ = 
2Snm(£'" + QmY^ \'\^m\d\Q\\'>i^n)'\ whcrc d%{g\+d\, Qx) = gxdX^, and we have used 
the eigenbasis gx = J2n 0n\ipn){ipn\- 

2. Gaussian states and a measure of non Gaussianity 

Let us consider a single-mode bosonic system described by the mode operator a 
with commutation relations [a, a^] = 1. A quantum state g is fully described by 
its characteristic function xiolW = Tr[£) Z)(A)] where D{X) = exp{Aa'l' — X*a} is the 
displacement operator. The canonical operators are given by g = (a + a^)/\/2 and 
p = (a — a^)/-\/2i with commutation relations given by [q,p] = i. Upon introducing the 
vector i? = (QtP)'^ , the covariance matrix cTg and the vector of mean values Xg of a 
quantum state g are defined as akj = ^{RkRj+RjRk)g~ {Rj) g{Rk) g and Xj = {Rj)g, 
where {0)g = Tr[£)0] is the expectation value of the operator O on the state g. A 
quantum state is said to be Gaussian if its characteristic function has a Gaussian 
form. Once the CM and the vectors of mean values are given, a Gaussian state is fully 
determined. 

The amount of nG 6[g] of a quantum state g may be quantified by the quantum 
relative entropy [34] S(^3|jTg) = Tr[f?(log p— log r)] between g and its reference Gaussian 
state Tg, which is a Gaussian state with the same covariance matrix cr as g. As for its 
classical counterpart, the KuUback-Leibler divergence, it can be demonstrated that < 
< 00 when it is definite, i.e. when suppp C suppr. In particular S(f3||Tg) = iff 
g = Tg [35l|36]. Since Tg is Gaussian Tr[(rg — p) logr^] = 0, i.e. S{g\\Tg) ~ S{Tg)~5{g) 
and we may write S[g] = S{Tg) — S{g) where S{g) = —TT[g \ogg] is the von Neumann 
entropy of g. Finally, since the von Neumann entropy of a single-mode Gaussian state 
may be written as h{y^dcta^) where h{x) = {x + ^) \og{x + i)) — {x — ^) log{x — ^) 
we have 

S[g] ^ hiy^d^g) - 5{g) . (1) 

A generic single-mode Gaussian state may be written as t ~ Us I'Ug where Us is 
a symplectic operation i.e. a unitary Us = exp —iH resulting from a Hamiltonian H at 
most quadratic in the field operators, and ;^ is a chaotic (maximum entropy) state with 
n„ = Tiliya^a] average thermal quanta, i.e. — ^kPk\^)(^\ ' Pk = '^^/(l + 'n^)^^^ in 
the Fock number basis. 
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3. Classical perturbations to a Gaussian state 

An infinitesimal perturbation of the eigenvalues pk of a Gaussian state r, i.e. 
Pk Pk + dpk results in a perturbed state g = J2kiPk + dpk)Us\k) {k\Ug which, 
in general, is no longer Gaussian. Since the nG of a state is invariant under 
symplectic operations we have S[q] = d[T]] where rj = U^qUs = '^^iPk + dpk)\k){k\ is 
diagonal in the Fock basis. The Gaussian reference of 77 is a thermal state with 
Tirj = n,y + dn = k {pk + dpk) average quanta and the nG may be evaluated upon 
expanding both terms in d[ri]{= S[g\) up to the second order, 



NonGaussianity of perturbed states is thus given by the sum of two contributions. 
The first term is the Fisher information of the probability distribution {pk}, which 
coincides with the classical part of the Bures distance in the Hilbert space. The 
second term is a negative contribution expressed in terms of the infinitesimal change 
of the average number of quanta. When traveling on surfaces at constant energy the 
amount of nG coincides with a proper distance in the Hilbert space and, in this case, 
it has a geometrical interpretation as the infinitesimal Bures distance. At the same 
time, since Bures distance is proportional to the QFI one, it expresses the statistical 
distinguishability of states, and we conclude that moving out from a Gaussian state 
towards its nouGaussian neighbours is a resource for estimation purposes. Similar 
conclusions can be made when comparing families of perturbations {dp} corresponding 
to the same infinitesimal change of energy dn^: in this case the different amounts 
of nG induced by the perturbations are quantified by the Bures distance minus a 
constant term depending on dn^ and the initial thermal energy rii/, i.e. the intial 
purity n = Trlg^] = Ti'[r'^] = {2ni, + We summarize the above statements in the 
following 

Theorem 1 // t\ is a Gaussian state and an infinitesimal variation of the value of A 
drives it into a state g\+dx with the same eigenvectors, then the QFI distance is equal 
to the nG S[g\+dx\ plus a term depending both on the infinitesimal variation of energy 
dn^ and on the initial purity 



In particular, for perturbations that leave the energy unperturbed the nG of the 
perturbed state coincides with the QFI distance, whereas, in general, it provides a 
lower bound. □ 

3.1. Examples of finite perturbations 

In order to explore specific directions in the neighbourhood of a Gaussian state r 
let us write the perturbation to the eigenvalues as dpk = e/ifc where {/ifc} is a given 
distribution. In this case the nG of the perturbed state is given by 




H{\)d\^ = 5[gx+dx] + 




{pk - iJ'kf 

2pk 



2nt{l + nt) 




) 



(3) 



\ fc 
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where /S.n^ — 'Y^^.^Pk^ ^k){k\a) a\k) Let us now consider the famihes of states generated 
by the convex combination g ~ {l — e)T + egf^ of the Gaussian states r with a the target 
state Qfj^ = ^kUs\k) {k\Ug , which itself is obtained by changing the eigenvahies pk of 
the initial Gaussian state to ^k- Again we exploit invariance of S[g\ under syniplectic 
operations and focus attention to the diagonal state 77 = U^gUs ~ (1 —e)rj + tri^ which 
has the same nG of g. This is the generalization to a finite perturbation of the analysis 
reported in the previous Section, and it is intended as a mean to find the maximally 
nonGaussian state obtainable starting from a given Gaussian. NonGaussianity of 
this kind of states can be written as 5[rf\ = S'(t^) — S{rii) = + 1/2] — Hlq^"^^], 
where H[q^^'>] = —^k'ik^^'^^l^k^ denotes the Shannon entropy of the distribution 
q^^^ = (1 — e)pk + e^fc and = Tr[rja'^a] = (1 — e)nt + en^ is the average number of 
quanta of 77. Notice that for a thermal state with rit quanta we have = h[nt + l/2]. 
By using the concavity of the Shannon entropy we obtain an upper bound for the nG 

S[v] < h[n,^ + 1/2] - h[nt + 1/2] + e(^h[nt + 1/2] - H[ii]j (4) 

In particular, if the two distributions have the same number of quanta, ri^ — nt, and 
thus nt = n^i, the bound on S[ri] only depends on the difference between the entropy 
of the initial and target distributions. 

Let us now consider perturbations towards some relevant distributions i.e. 

Poissonian /i^'''' = lif'^""'^'; thermal /i^*'', and Fock /i^,'^'' ~ Sk,n^, and evaluate the 
nG of states obtained as convex combination of a thermal state with nt quanta and 
a diagonal quantum state with a Poissonian, thermal or Fock distributions and 
quanta. 




Figure 1. NonGaussianity S[ti] as a function of the energy of the target state 
rip , with = 4 and for different values of e = 0.3, 0.7, 0.9 (red, green and blue, 
respectively) for Poissonian target (upper loft panel), thermal (upper right) and 
Fock (lower left). The lower right panel shows the nG Slrj] as a function of e, 
with Ufj, = nt = 4 and for convex combination with a Poissonian (red) and a Fock 
(blue) state. Dashed lines correspond to the expansion of nG at the second order 
in e, as obtained in Eq. l[3} for the corresponding perturbations. 



In Fig. [T] we plot the nG of the convex combination g as a function of for 
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different values of e: If we consider the eonvex combination with a Fock state, the nG 
simply increases monotonically with the energy of the added state. For combinations 
with Poissonian and thermal distributions we have a maximum for — > 0, then a 
local minimum (which for the thermal distribution corresponds trivially to 5[q\ = 0) 
for — nt, and then the nonG increases again for higher values of nx . This implies 
that in order to increase nG the best thing to do is to perturb the initial Gaussian 
either with a highly excited state ^ 1 or with the vacuum state |0) (0|. If wc choose 
Uq — Uf^ = ut wc know from Eq. ([2]) that nG for small perturbations is equivalent to 
the Bures infinitesimal distance between the two probability distributions fik and pk ■ 
In Fig. [TJi we show the nG as a function of e: As it is apparent from the plot the nG 
obtained by adding a Fock state is always much larger that the one for a Poissonian 
profile. We observe that in the latter case the expansion at the second order obtained 
in Eq. is still accurate for values of e approaching 1, while it fails to be accurate 
for e > 0.2 for a Fock state, becoming an upper bound on the exact amount of nG. We 
have also investigated what happens by considering a target distribution randomly 
chosen on a finite subspace of the infinite Hilbert space: Again we have obtained that, 
at fixed energy of the target state, perturbing with a Fock state yields the biggest 
increase of nG. This may be easily generalized to a system of d bosonic modes, where 
the most general Gaussian state is described by d{2d + 3) independent parameters. 

4. Perturbations at fixed covariance matrix 

Gaussian states are known to be extremal state at fixed covariance matrix for several 
relevant quantities, e.g. channel capacities and entanglement measures [37] . Therefore, 
one may wonder whether perturbing a Gaussian state at fixed covariance matrix 
may be quantified in convenient way for the purposes of quantum estimation. This 
indeed the case: the nG provides an upper bound to the QFI distance at fixed 
covariance matrix and thus have an operational interpretation in terms of statistical 
distinguishability. This is more precisely expressed by the following theorem [1] 

Theorem 2 // tx is a Gaussian state and an infinitesimal variation of the value of 
A drives it into a state gx+d\ with the same covariance matrix, then the nG S[gx+dx\ 
provides an upper bound to the QFI distance 

H{X)d\^ < 6[gx+dx] ■ 

Proof: If gx+d\ and tx have the same CM then the nG of gx+d\, S[gx+d\] = 
S(gA+dA I ^a) = (A)dA^ , where the so-called Kubo-Mori-Bogolubov information H{X) 
|381 139] provides an upper bound for the quantum Fisher information H{X) < H{X) 
[40] , thus proving the theorem. □ 

The above theorem says that a larger nG of the perturbed state may correspond 
to a greater distinguishability from the original one, thus allowing a more precise 
estimation. Of course, this is not ensured by the theorem, which only provides an 
upper bound to the QFI. One may wonder that when gx+d\ is itself a Gaussian state 
the theorem requires H{X) = 0, i.e. no reliable estimation is possible. Indeed, this 
should be the case, since Gaussian states are uniquely determined by the first two 
moments and thus the requirement that the perturbed ta+^a and the original state ta 
are both Gaussian and have the same covariance matrix implies that they are actually 
the same quantum state. 
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5. Conclusions 

In conclusion, we have addressed the nG of states obtained by weakly perturbing a 
Gaussian states and have investigated the relationships with quantum estimation. We 
found that nG provides a lower bound to the QFI distance for classical perturbations, 
i.e. perturbations to eigenvalues leaving the eigenvectors unperturbed, and an upper 
bound for perturbations leaving the covariance matrix unperturbed. For situations 
where the CM is changed by the perturbation we have no general results. On the other 
hand, it has been already shown that non-Gaussian states improve quantum estimation 
of both unitary perturbations as the displacement and the squeezing parameters [JT] 
and nonunitary ones as the loss parameter of a dissipative channel [42j . Overall, 
our results show that the geometry of nonGaussian states in the neighbourhood of 
a Gaussian state is definitely not trivial and cannot be subsumed by a differential 
structure. Despite this fact, the analysis of perturbations to a Gaussian state may 
help in revealing when, and to which extent, nG is a resource for quantum estimation. 
We have also analyzed the nG of specific families of perturbed Gaussian states with 
the aim of finding the maximally non Gaussian state obtainable from a given Gaussian 
one. 
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